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Abstract 



We study the power spectrum of the mass density perturbations in an 
inflation scenario that includes thermal dissipation. We show that the con- 
dition on which the thermal fluctuations dominate the primordial density 
perturbations can easily be realized even for weak dissipation, i.e., the rate of 
dissipation is less than the Hubble expansion. We find that our spectrum of 
primordial density perturbations follows a power law behavior, and exhibits 
a "thermodynamical" feature - the amplitude and power index of the spec- 
trum depend mainly on the thermodynamical variable M, the inflation energy 
scale. Comparing this result with the observed temperature fluctuations of 
the cosmic microwave background, we find that both amplitude and index of 

the power spectrum can be fairly well fitted if M ~ 10^^ — 10^^ GeV. 
PACS number(s): 98.80.Cq, 98.80.Bp, 98.70.Vc 
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I. INTRODUCTION 



In the past decade, there has been a number of studies on dissipative processes associated 
with the inflaton decay during its evolution. These studies have shed hght into the possible 
effects of the dissipative processes. For instance, it was realized that dissipation effectively 
slows down the rolling of the inflaton scalar fleld (f) toward the true vacuum. These processes 
are capable of supporting the scenario of inflation ||T],0|. Recently, inspired by several new 
developments, the problem of inflation with thermal dissipation has attracted many re- 
investigations. 

The flrst progress is from the study of the non-equihbrium statistics of quantum flelds, 
which has found that, under certain conditions, it seems to be reasonable to introduce a 
dissipative term (such as a friction-like term) into the equation of motion of the scalar fleld 
(f) to describe the effect of heat contact between the fleld and a thermal bath. These 
studies shown that the thermal dissipation and fluctuation will most likely appear during 
the inflation if the inflaton is coupled to light flelds . However, to realize sufficient e-folds 
of inflation with thermal dissipation, this theory needs to introduce tens of thousands of 
scalar and fermion flelds interacting with the inflaton in an ad hoc manner Namely, 
it is still far from a realistic model. Nevertheless, this study indicates that the condition 
necessary for the "standard" reheating evolution - a coupling of inflaton with light flelds - 
is actually also the condition under which the effects of thermal dissipation during inflation 
should be considered. 

Secondly, in the case of a thermal bath with a temperature higher than the Hawking tem- 
perature, the thermal fluctuations of the scalar fleld plays an important and even dominant 
role in producing the primordial perturbations of the universe. Based on these results, the 
warm inflation scenario has been proposed. In this model, the inflation epoch can smoothly 
evolve to a radiation-dominated epoch, without the need of a reheating stage Dynam- 
ical analysis of systems of inflaton with thermal dissipation gives further support to this 
model. It is found that the warm inflation solution is very common. A rate of dissipation 
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as small as 10~^ if, H being the Hubble parameter during inflation, can lead to a smooth 
exit from inflation to radiation. 

Warm inflation also provides explanation to the super-Hubble suppression. The standard 
inflationary cosmology, which is characterized by an isentropic de Sitter expansion, predicts 
that the particle horizon should be much larger than the present-day Hubble radius c/Hq. 
However, a spectral analysis of the COBE-DMR 4-year sky maps seems to show a lack of 
power in the spectrum of the primordial density perturbations on scales equal to or larger 
than the Hubble radius c/ Hq [I],^. A possible explanation of this super-Hubble suppression is 
given by hybrid models, where the primordial density perturbations are not purely adiabatic, 
but mixed with an isocurvature component. The warm inflation is one of the mechanisms 
which can naturally produce both adiabatic and isocurvature initial perturbations 0. 

In this paper, we study the power spectrum of mass density perturbations caused by 
inflation with thermal dissipation. One purpose of developing the model of warm inflation 
is to explain the amplitudes of the initial perturbations. Usually, the amplitude of initial 
perturbations from quantum fluctuation of the inflaton depends on some unknown param- 
eters of the inflation potential. However, for the warm inflation model, the amplitude of 
the initial perturbations is found to be mainly determined by the energy scale of inflation, 
M. If M is taken to be about 10^^ GeV, the possible amplitudes of the initial perturbations 
are found to be in a range consistent with the observations of the temperature fluctuations 
of the cosmic microwave background (CMB) ]10|]. That is, the thermally originated initial 
perturbations apparently do not directly depend on the details of the inflation potential, 
but only on some thermodynamical variables, such as the energy scale M. This result is not 
unexpected, because like many thermodynamical systems, the thermal properties including 
density fluctuations should be determined by the thermodynamical conditions, regardless of 
other details. 

Obviously, it would be interesting to find more "thermodynamical" features which con- 
tain only observable quantities and thermodynamical parameters, as these predictions would 
be more useful for confronting models with observations. Guided by these considerations, 
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we will extend the above-mentioned qualitative estimation of the order of the density per- 
turbations to a quantitative calculation of the power spectrum of the density perturbations. 
We show that the power spectrum of the warm inflation does not depend on unknown pa- 
rameters of the inflaton potential and the dissipation, but only on the energy scale M. The 
spectrum is found to be of power law, and the index of the power law can be larger or less 
than 1. More interestingly, we find that for a given M, the amplitude and the index of the 
power law are not independent from each other. In other words, the amplitude of the power 
spectrum is completely determined by the power index and the number M. Comparing this 
result with the observed temperature fluctuations of the CMB, we find that both amplitude 
and index of the power spectrum can be fairly well fitted if M ~ 10^^ — 10^^ GeV. 

This paper is organized as follow: In Sec. II we discuss the evolution of the radiation 
component for inflationary models with dissipation prescribed by a field-dependent friction 
term. In particular, we scrutinize the physical conditions on which the thermal fluctuations 
dominate the primordial density perturbations. Section III carries out the calculations of 
the power spectrum of the density perturbations of the warm inflations. And finally, in Sec. 
IV we give the conclusions and discuss further observational tests. 



II. INFLATION WITH THERMAL DISSIPATION 
A. Basic equations 

Let us consider a flat universe consisting of a scalar inflaton field and a thermal bath. 
Its dynamics is described by the following equations p. The equations of the expanding 
universe are 
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where H = R/R is the Hubble parameter, and mpi = the Planck mass. V{(j)) is the 

effective potential for field 0, and pr is the energy density of the thermal bath. Actually the 
scalar field (f) is not uniform due to fiuctuations. Therefore, the field in Eqs.(2.1) and (2.2) 
should be considered as an average over the fiuctuations. 

The equation of motion for scalar field in a de-Sitter universe is 

+ 3i70 + 10 - e-^^*VV + V\(t)) = 0, (2.3) 

where the friction term r<p describes the interaction between the field and a heat bath. 
Obviously, for a uniformed field, or averag ed 0, the term of Eq. (2.3) can be ignored. 
Statistical mechanics of quantum open systems has shown that the interaction of quantum 
fields with thermal or quantum bath can be described by a general fiuctuation-dissipation 
relation It is probably reasonable to describe the interaction between the infiaton and 



the heat bath as a "decay" of the infiaton [|T2|. These results support the idea of introducing 



a damping or friction term into the field equation of motion. In particular, the friction term 
with the form in Eq. (2.3), r0, is a possible approximation for the dissipation of field in 
a heat bath environment in the near-equilibrium circumstances. In principle, F can be a 
function of 0. In the cases of polynomial interactions between field and bath environment, 
one may take the polynomial of for F, i.e., F = Fm0™. The friction coefficient must be 
positive definite, hence F^ > 0, and the dissipative index of friction m should be zero or 
even integer if V{(j)) is invariant under the transformation — > —0. 

The equation of the radiation component (thermal bath) is given by the first law of 
thermodynamics as 

Pr + AHpr = F0^ (2.4) 

The temperature of the thermal bath can be calculated by p,. = {n^ /SO) gesT'^, ges being the 
effective number of degrees of freedom at temperature T. 

The warm infiation scenario is generally defined by a characteristic that the thermal 
fiuctuations of the scalar field dominate over the quantum origin of the initial density per- 
turbations. Because the thermal and quantum fiuctuations of the scalar field are proportional 



to T and H respectively, a necessary condition for warm inflation models is the existence of 
a radiation component with temperature 



during the inflationary expansion. Eq. (2.5) is also necessary for maintaining the thermal 
equilibrium of the radiation component. In general, the time scale for the relaxation of a 
radiation bath is shorter for higher temperature. Accordingly, to have a relaxing time of the 
bath shorter than the expansion of the universe, a temperature higher than H is generally 
needed. 

As a consequence of Eq. (2.5), warm inflation scenario requires that the solutions of Eqs. 
(2.1) - (2.4) should contain an inflation era, followed by smooth transition to a radiation- 
dominated era. Dynamical system analysis also confirmed that for a massive scalar field 
Vi^cj)) = |M^0^, the warm inflation solution of Eqs. (2.1) - (2.4) is very common. A smooth 
exit from inflation to radiation era can be established even for a dissipation with F as small 
as 10^^ H 1^. A typical solution of warm inflation will be given in next section. 



Since warm inflation solution does not rely on a specific potential, we will employ the 
popular 0^ potential commonly used for the "new" inflation models. It is 



To have slow-roll solutions, the potential should be flat enough, i.e., A < (M/mpi)^, where 



For models based on the potential of Eq. (2.6), the existence of a thermal component 
during inflation seems to be inevitable. In order to maintain the field close to its minimum 
at the onset of the inflation phase transition, a thermal force is generically necessary. In 
other words, there is, at least, a weak coupling between field and other fields contributing 



T > H 



(2.5) 



B. Evolution of Radiation component during inflation 




(2.6) 



V{0) = = Xa\ 
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to the thermal bath. During the slow roll period of inflation, the potential energy of the (f) 
field is fairly constant, and their kinetic energy is small, so that the interaction between the 
4> field with the fields of the thermal bath remains about the same as at the beginning. As 
such, there is no compelling reason to ignore these interactions. 

Strictly speaking, we should use a finite temperature effective potential 1^(0, T). How- 
ever, the correction due to finite temperature is negligible. The leading temperature correc- 
tion of the potential (2.6) is AT^^^. On the other hand, as mentioned above, we have A < 
(M/mpi)^ for the flatness of the potential. Therefore, AT^ < MVm|,i ~ {M/mpifH^ < H^, 
i.e., the influence of the flnite temperature effective potential can be ignored when < mpi. 

Now, we try to flnd warm inflation solutions of Eqs. (2.1) - (2.4) for weak friction T < H. 
In this case, Eqs. (2.1) - (2.3) are actually the same as the "standard" new inflation model 
when 

Pr « V^(O). (2.7) 
Namely, we have the slow-roll solution as 

^ n^L^m .28) 

« V{0), (2.9) 

and 



where the subscript i denotes the starting time of the inflation epoch. 

During the stage of -C cr, it is reasonable to neglect the 0^ term in Eq. (2.3). We have 
then 



4)+{3H + r)0 - 4Xa^(j) = 0. 

Considering F < H, an approximate solution of can immediately be found as 
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(2.11) 



(f) = (f)ie' 



(2.12) 



where a ~ 



A^/^(mpi/M)^/27r and 0i is the initial value of the scalar field. 



Substituting solution (2.12) into Eq. (2.4), we have the general solution of (2.4) as 



where A = Q;^ifr^0^+^/[(m + 2)a + 4], S = ^^(0) - A, and Pr(0) is the initial radiation 
density. Obviously, the term B in Eq. (2.12) describes the blowing away of the initial 
radiation by the inflationary exponential expansion, and the term A is due to the generation 
of radiation by the field decay. 

According to Eq. (2.13), the evolution of the radiation has two phases. Phase 1 covers 
the period during which the B term is dominant, and radiation density drops drastically 
due to the infiationary expansion. The component of radiation evolves into phase 2 when 
the A term becomes dominant, where the radiation density increases due to the friction of 
the (p field. Namely, both heating and infiation are simultaneously underway in phase 2. 
Therefore, this phase is actually the era of infiation plus reheating. 

The transition from phase 1 to phase 2 occurs at time tb determined by {dp/dt)t^ = 0. 
We have 



-4Ht 



(2.13) 




(2.14) 



where a = (7r^/30)gfeff- Then the radiation density at the rebound time becomes 




(2.15) 



From Eqs. (2.12) and (2.13), the radiation density in phase 2 is given by 




(2.16) 



Since H ~ (M/mpi)M, Eq. (2.16) can be rewritten as 




(2.17) 
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On the other hand, from (2.12), we have 

>-AV^(^)^(M)V(0). (2.18) 

Therefore, in the case of weak dissipation T < H, we have 

Pr{t) < 4>{tf/2. (2.19) 

This is consistent with the condition of inflation Eq. (2.7) when Eq. (2.9) holds. 

Eqs. (2.7) and (2.9) indicate that the inflation will come to an end at time tf when the 
energy density of the radiation components, or the kinetic energy of field, 0^/2, become 
large enough, and comparable to V{Qi). From Eqs. (2.17) and (2.18), tf is given by 

In general, at the time when the phase 2 ends, or a radiation-dominated era starts, the 
potential energy may not be fully exhausted yet. In this case, a non-zero potential V will 
remain in the radiation-dominated era, and the process of decaying into light particles is 
still continuing. 

However, considering A^/^(mpi/M)^(r/i7) < 1, the right hand side of Eq. (2.17) will 
always be less than 1 when 0(i) is less than a. This means that, for weak dissipation, phase 
2 cannot terminate at (j){t) < a, or V^(0(t/)) 7^ 0. Therefore, under weak dissipation, phase 
2 will end at the time tf when the potential energy V{4>) is completely exhausted, i.e., 

(p{tf) - a. (2.21) 

This means that no non-zero V remains once the inflation exits to a radiation-dominated 
era, and the heating of (p decay also ends at tf. 



C. Temperature of radiation 



From Eq. (2.13), one can find the temperature T of the radiation in phases 1 {t < tb) 
and 2 {t > tb) as 
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, , , T,e-^(*-*^), if t<h, ^ ^ 

T{t) = { (2.22) 

rp^^{m+2)aH{t-t,)/A^ if tf > t > U, 



where 



The temperature Tf at the end of phase 2 is 



m + 2 
4 



(2.23) 



= T(t/) = Tbe^rn+2)aH{tf-t,)/4^ ^2.24) 

where is given by Eq. (2.21). 

Since T{t) is increasing with t in phase 2, the condition (2.5) for warm inflation can be 
satisfied if T(t/) > H, or 

Pritf) > aHf. (2.25) 
Using Eqs. (2.17) and (2.21), condition (2.25) is reahzed if 

r ^ / (T \2 / M y 

H Vmpi/ Vmpi/ 

Namely, T can be as small as IQ-^'^H for M ~ 10^*^ GeV, and a ~ 10^^ GeV. Therefore, 
the radiation solution (2.13), or warm inflation, should be taken into account in a very wide 
range of dissipation 

IQ-^^H <T <H. (2.27) 

This result is about the same as that given by dynamical system analysis : a tiny friction 
r may lead the infiaton to a smooth exit directly at the end of the inflation era. 

A typical solution of the evolution of radiation temperature T(t) is demonstrated in Fig. 
1, for which parameters are taken to be M = 10^^ GeV, a = 2.24 ■ 10^^ GeV, T2 = lO'^Hi 
and Qcs = 100. Actually, (^off-factor is a function of T in general. However, as can be 
seen below, the unknown function gea{T) has only a slight effect on the problems under 
investigation. Figure 1 shows that the rebound temperature can be less than H. In this 
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case, the evolution of T{t) in phase 2 can be divided into two sectors: T < H ior t < te, and 
T > H ior t > te, where tg is defined by T{te) = H. We should not consider the solution 
of radiation to be physical iiT < H since it is impossible to maintain a thermalized heat 
bath with the radiation temperature less than the Hawking temperature H of an expanding 
universe. Nevertheless, the solution (2.13) should be available \lt> tg. Therefore, one can 
only consider the period of te < ^ < ^/ as the epoch of the warm inflation. 

Figure 1 also plots the Hubble parameter H{t). The evolution of H{t) is about the same 
as in the standard new inflation model, i.e., H(t) ~ Hi in both phases 1 and 2. In Fig. 1, it 
is evident that the inflation smoothly exits to a radiation era at t/. The Hubble parameter 
H{t) also evolves from the inflation H{t) ~ constant to a radiation regime H{t) oc t~^. 

The duration of the warm inflation is represented by (t/ — te) then. The number of 
e-folding growth of the comoving scale factor R during the warm inflation is given by 

N= Hdt^- — In^. (2.28) 

Jte {m + 2)a H ^ ^ 

One can also formally calculate the number of e-folds of the growth in phase 2 as 

N2 = Hdt^ — In ^, (2.29) 

4 {m + 2)a Tb' ^ ^ 

and the number of e-folds of the total growth as 

r*/ „ , 4 , /T- 



N, 



r Hdt^ In (^) + Ht,. (2.30) 

Jo (m + 2)a \TbJ 



It can be found from Eqs. (2.28) - (2.30) that both N2 and Nt depend on the initial 
value of the fleld 0^ via T^, but N does not. The behavior of T at the period t > te is 
completely determined by the competition between the diluting and producing radiation 
at t > tb. Initial information about the radiation has been washed out by the inflationary 
expansion. Hence, the initial (pi will not lead to uncertainty in our analysis if we are only 
concerned the problems of warm evolution at the period te < t < tf. 

III. THE PRIMORDIAL DENSITY PERTURBATIONS 
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A. Density fluctuations of the 4> field 



The fluctuations of field can be calculated by the similar way as stochastic inflations 



Recall that the coarse-grained scalar field is actually determined from the decompo- 



sition between background and high frequency modes, i.e. 

<l>(x,t) =0(x,t) + g(x,t), (3.1) 

where $(x, t) is the scalar field satisfying 

$ + 3i7$ - e~^^*V^$ + V'($) = 0. (3.2) 

g(x, t) in Eq.(3.1) contains all high frequency modes and gives rise to the thermal fluctua- 
tions. Since the mass of the field can be ignored for the high frequency modes, we have 

g(x, t) = j (fkWiM) [akak(t)e-'''-" + 4cr*(t)e''^-"] (3.3) 

where k is comoving wave vector, and modes crk(t) is given by 

e-'^\ (3.4) 



1 1 

t^k(t) 



H 

(27r)3/2 L^^~'"fc, 

and r = —H~^ exp{—Ht) is the conformal time. Eq.(3.3) is appropriate in the sense that the 
self-coupling of the field is negligible. Considering the high frequency modes are mainly 
determined by the heat bath, this approximation is reasonable. The window function iy(|k|) 
is properly chosen to filter out the modes at scales larger than the horizon size H~^, i.e., 
W{k) = 6{k — kh(t)), where kh{t) ~ (l/7r)if exp(ift) Qis the lower limit to the wavenumber 
of thermal fluctuations. 

From Eqs.(3.1) and (3.3), with the slow-roll condition, Eq.(3.2) renders 



^The coefficient l/vr actually depends on the details of the cut-off function, which may not be 
step-function-like. For instance, considering causality, the cut-off function can be soft, and the 
longest wavelength of fluctuations can be a few times of the size of horizon Q 
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(3.5) 



and 



( 




g(x,t). 



(3.6) 



Eq.(3.5) can be rewritten as 



d(j){'x, t) 
Jt 



1 5F[<P{^,t)] 
3H (50 



+ ^(x,t) 



(3.7) 



where 



F[0] = Jd'^ [-(e-^*V0)2 + T/(0) 



(3.8) 



Eq. (3.7) is, in fact, the rate equation of the order parameter of a system with free energy 
F[(j)]. It describes the approach to equihbrium for the system during phase transition. 

Using the expression of free energy (3.8), the slow-roll solution (2.8) can be rewritten as 



Hence, in the case of weak dissipation (F < H), Eq. (3.7) is essentially the same as the 
slow-roll solution (2.8) or Eq. (3.9) but with fluctuations rj. The existence of the noise 
field ensures that the dynamical system properly approaches the global minimum of the 
inflaton potential V{(f)). Strictly speaking, both the dissipation F and fluctuations rj are 
consequences derived from g(x, t). They should be considered together. However, it seems 
to be reasonable to calculate the fluctuations alone if the dissipation is weak. 

Unlike (3.2), the Langevin equation (3.7) is of first order (0) due to the slow-roll con- 
dition. Generally, thermal fluctuations will cause both growing and decaying modes [Q] 0. 
Therefore, the slow-roll condition simplifies the problem from two types of fluctuation modes 
to one, i.e., we can directly calculate the total fluctuation as the superposition of various 



^We thank the referee for pointing this problem out. 



d(j) 



1 dF[(p] 
3H + T d0 



(3.9) 
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fluctuations. It has been shown [|T^] that during the eras of dissipations, the growth of the 
structures in the universe is substantially the same as surface roughening due to stochastic 
noise. The evolution of the noise-induced surface roughening is described by the so-called 
KPZ-equation [jT6[. Eqs.(3.5) or (3.7), which includes terms of non-linear drift plus stochastic 
fluctuations, is a typical KPZ-like equation. 

From Eq.(3.6), the two-point correlation function of f]{'x,t) can be found as 



(r/(x,t)r7(x',t')) 



47r2 



1 + 



sin(fcfc|x - x'l) 
kh |x — x'l 



6{t-t'), 



(3.10) 



exp{H/TTT) - 1 

where l/[exp(if/7rT) — 1] is the Bose factor at temperature T. Therefore, when T > H, we 
have 



(r/(x,t)77(x,t')) 



5{t-t'). 



(3.11) 



This result can also be directly obtained via the fluctuation-dissipation theorem ||T7| , p!8| . In 
order to accord with the dissipation terms of Eq. (3.7), the fluctuation-dissipation theorem 
requires the ensemble average of r] to be given by 



(r/) = 



(3.12) 



and 



(r7(x,t)r/(x,t')) = ^5(t-0 



(3.13) 



The variance D is determined by 



D = 2- 



T 



(3.14) 



where U = {A7r/3)H~^ is the volume with Hubble radius H^^. In the case of weak dissipa- 
tion, we then recover the same result as in Eq.(3.11), 



D = H'^T/2n. 



(3.15) 



When T = H, we obtain 
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D = -, (3.16) 
which agrees exactly the result derived from quantum fluctuations of (p-field . Therefore, 



the quantum fluctuations of inflationary (j) field are equivalent to the thermal noises stim- 
ulated by a thermal bath with the Hawking temperature H. Eqs. (3.15) and (3.16) show 
that the condition (2.5) is necessary and sufficient for a warm inflation. 

For long- wavelength modes, the V'{(j)) term is not negligible. It may lead to a suppression 
of correlations on scales larger than |y (0)1"^/^. However, before the inflaton actually rolls 
down to the global minimum, we have |^"(0 ^ o")|^^/^ > H^^. The so-called abnormal 



dissipation of density perturbations [0 may produce more longer correlation time than 
H. Therefore in phase 2, i.e., the warm inflation phase H < T < M, the long-wavelength 
suppression will not substantially change the scenario presented above. 

The fluctuations 6(f) of the field can be found from linearizing Eq. (3.7). If we only 
consider the fluctuations 6(j) crossing outside the horizon, i.e., with wavelength ~ H^^, the 
equation of 6(j) is 

-dr = - 3H+T ^'^+^- 

For the slow- roll evolution, we have |y"(0)| <^ 9H^ 0. One can ignore the V" {(j)) term on 
the right hand side of Eq. (3.17). Accordingly, the correlation function of the fluctuations 
is 

(Smmt')) - D^^e-(*-*')^^/(™, t > t', (3.18) 

hence 

~ ^HT (3.19) 

Thus, in the period <t <tf the density perturbations on large scales are produced by the 
thermal fluctuations that leave the horizon with a Gaussian-distributed amplitude having a 
root-mean-square dispersion given by Eq. (3.19). 
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Principally, the problem of horizon crossing of thermal fluctuations given by Eq. (3.7) is 
different from the case of quantum fluctuations, because the equations of H and iJ, (2.1) and 
(2.2) contain terms in p^- However, these terms are insignificant for weak dissipation [Eq. 
(2.19)] in phase 2. Thus Eqs.(2.1) and (2.2) depend only nominally on the evolution of pr- 
Accordingly, for weak dissipation, the behavior of thermal fluctuations at horizon crossing 
can be treated by the same way as the evolutions of quantum fluctuations in stochastic infla- 
tion. In that theory, quantum fluctuations of inflaton are assumed to become classical upon 
horizon crossing and act as stochastic forces. Obviously, this assumption is not necessary for 
thermal fluctuations. Moreover, we will show that in phase 2 the thermal stochastic force 
HT is contingent upon the comoving scale of perturbations by a power law [Eqs. (2.21) and 
(3.21)], and therefore the power spectrum of the thermal fluctuations obeys the power law. 
This make it more easier to estimate the constraint quantity in the super-horizon regime. 

Accordingly, the density perturbations at the horizon re-entry epoch are characterized 

by i 

C-l] ^^^n^. (3.20) 
VpA 02 + (4/3)p. 

All quantities in the right-hand side of Eq. (3.20) are calculated at the time when the 

relevant perturbations cut across the horizon at the inflationary epoch. 

Using the solutions of and pr of warm inflation (2.12) and (2.13), Eq. (3.20) gives 

(5p\ / 5.33^/2+4 /^^\^ /T\Ki7Tl) 

l^p " 1^2-+3.7rW2+3j \g,Qa--) \h) ' ^ ^ 

where the dimensionless parameter 7^ = r^-ff™"^, and T is the temperature at the time 

when the considered perturbations 5pr crossing out of the horizon R-^ ~ Hr\ Eq. (3.21) 

shows that the density perturbations are insensitive to the (yfcg-factor. 



B. Power law index 



Since inflation is immediately followed by the radiation dominated epoch, the comoving 
scale of a perturbation with crossing over (the Hubble radius) at time t is given by 
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where Tq and Hq are the present CMB temperature and Hubble constant respectively. Eq. 
(3.22) shows that the smaller t is, the smaller k will be. This is the so-called "first out - last 
in" of the evolution of density perturbations produced by the inflation. 

Using Eqs. (2.22) and (3.22), the perturbations (3.21) can be rewritten as 

(^(j-J^ oc)t("*-6W^ if k>ke, (3.23) 

where kg is the wavenumber of perturbations crossing out of horizon at te- It is 

k, = 27ri/^e^(*-*/) ~ 27r//^e-^. (3.24) 

Therefore, the primordial density perturbations produced during warm inflation are of power 
law with an index (m — 6)a/4. We may also express the power spectrum of the density 
perturbations at a given time t. It is 



^ ocA;2+", if k>ke, (3.25) 



n^l+ ('^^) a. (3.26) 



where the spectral index n is 



4 

Clearly, for m — 6, the warm inflation model generates a flat power spectrum n — 1, yet 
the power spectrums will be tilted for m ^ Q. The dissipation models F = Tm4'"^ may not 
be realistic for higher m, but we will treat m like a free parameter in order to show that the 
results we concerned actually are not very sensitive to these parameters. 

The warm inflation scenario requires that all perturbations on comoving scales equal to 
or less than the present Hubble radius originate in the period of warm inflation. Hence, the 
longest wavelength of the perturbations (3.24), i.e., 27r/ke, should be larger than the present 
Hubble radius Hq^. We have then 
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where we have used {Tq/Hq) > (Tf/H), asTf <M. Using Eq. (2.28), the condition (3.27) 
gives an upper bound to a for a given m as 

4 N HTf/H) 



(3.28) 



n 



(3.29) 



.m + 2j ln{To/Ho)' 
Thus, the possible area of the index n can be found from Eq. (3.27) as 

1 — (6 — m)Q;niax/4 to 1, if m < 6, 
1 to 1 + (m — 6)Q;niax/4, if m > 6. 

Therefore, the power spectrum is positive-titled (i.e., n > 1) if m > 6, and negative-titled 
(n < 1) if m < 6. Figure 2 plots the allowed area of n as a function of the inflation mass 
scale M. Apparently, for M > 10^^ GeV, the tilt |n — 1| should not be larger than about 
0.15 regardless of the values of m from 2 to 12. 



C. Amplitudes of perturbations 



To calculate the amplitude of the perturbations we rewrite spectrum (3.25) into 

2\ / , \ n-1 

, if k>ke, (3.30) 



where ko — 2t:Hq. A is the spectrum amplitude normalized on scale k — kg, corresponding 
to the scale on which the perturbations re-enter the Hubble radius I/Hq at present time. 
From Eqs. (3.21), and (3.23), we have 

' 5.33m/2+4 \^/ /HoTfV^ /T\^ 



A 



^{n-l)H{tf-t) 



(3.31) 



Applying Eq. (2.21), the radiation temperature at the moment of horizon-crossing, t, can 
be expressed as T{t) = Tf exp[(m + 2)aH{t-t /) /4] . With the help of Eq. (2.28) , we obtain 



1 + I : — a 



liT' (I) " <='"-'"""-" = exp I 



On the other hand, using Eqs. (2.20), (2.23) and (2.28), one has 



N 



(3.32) 



In 



30 Vmpi 



Q m 

a ^ ^ . 



(3.33) 



Substituting Eqs. (3.32) and (3.33) into Eq. (3.31), we have finally 

•*=(64?Tf) i;^) ii) a-»exp{(„-l)[l+(— )„]iv}. (3.34) 

Eq. (3.34) shows that the amplitude A does not contain the unknown gfeff-factor. More- 
over, a can be expressed by n and m through Eq. (3.26), and can be expressed by a and 
M via Eq. (2.28). Therefore, the amplitude of the initial density perturbations. A, is only 
a function of M, n, and m. 

Figures 3 and 4 plot the relations between the amplitude A and index n for various 
parameters M and m. In the case of m = 6, n = 1, the relation of A and a is plotted in 
Fig. 5. It can be seen from Figs. 3, 4 and 5 that for either m > 6 or m < 6, the amplitude 
A is significantly dependent on M, but not so sensitive to m. Namely, the testable A-n 
relationship is mainly determined by a thermodynamical variable, the energy scale M. This 
is a "thermodynamical" feature. The relationship between A and plotted in Figs. 6 and 
7 also show this kind of "thermodynamical" feature: the A-N relation depends mainly on 
M. 

For comparison, the observed results of A and n derived from the 4-year COBE-DMR 
data (quadrupole moment Qrms-ps ~ I^.^kJ^Ii^K ^"^^ n ~ 1.2 ± 0.3 [|lOl) are plotted in 
Figs. 3, 4 and 5. The observationally allowed A-n range is generally in a good agreement 
with the predicted A-n curve if M ~ 10^^ — 10^^ GeV, regardless the parameter m. Figures 
3 and 4 also indicate that if the tilt of spectrum |n — 1| is larger than 0.1, the parameter 
area of M < 10^^ GeV will be ruled out. Therefore, the warm inflation seems to fairly well 
reconcile the initial perturbations with the energy scale of the inflation. 

IV. CONCLUSIONS AND DISCUSSION 

Assuming that the infiaton </)-field undergoes a dissipative process with F^^, we have 
studied the power spectrum of the mass density perturbations. In this analysis, we have 
employed the popular 0^ potential. However, only one parameter, the mass scale of the in- 
flation M , is found to be important in predicting the observable features of power spectrum, 
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i.e., the amplitude A and index n. Actually, the warm inflation scenario is based on two 
thermodynamical requirements: (a) the existence of a thermalized heat bath during infla- 
tion, and (b) that the initial fluctuations are given by the fluctuation-dissipation theorem. 
Therefore, we believe that the "thermodynamical" features - A and n depend only on M 
- would be generic for the warm inflation. This feature is useful for model testing. Hence, 
the warm inflation can be employed as an effective working model when more precise data 
about the observable quantities A, n etc. become available. The current observed data of 
A and n from CMB are consistent with the warm inflation scenario if the mass scale M of 
the inflation is in the range of 10^^ — 10^^ GeV. 
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FIGURES 

FIG. 1. A typical solutions of the evolutions of ^ field and radiations in warm inflation in 

which r = r2</>^, V{(t>) = \{4P' — cr'^)'^ and = M^, M being the inflaton energy scale. The 
parameters are taken to be M = lO^^GeV, a = 2.24 • 10^^ GeV, = IQ'^Hi and c/eff = 100. The 
dot-dashed and solid lines are for H{t) and the radiation temperature T respectively. 4 is the time 
at which the temperature rebound (T;,) and t(. the time of T = H. The inflation ends at when 
the temperature is Tf. T and H are in units of Hi = [87rF(0)/3mpj]^/^, and t is in units of H^^. 

FIG. 2. The allowed area of power law index n as a function of the mass scale M for various 
m. For given M and m, the possible n should lie between the line n = 1 and the corresponding 
curve of m. For m = 6, the only solution is n = 1. 

FIG. 3. The amplitudes of the power spectrum as a function of n in the area of n < 1. The 
mass scales M are labeled at the curves. The dotted, dashed, and solid lines are for m = 0, 2 and 
4 respectively. All curves end at the points when the corresponding warm inflation durations N 
are less than 55. The region within the dot-dashed box is the allowed area of (n. A) given by the 
4-year COBE-DMR data. 

FIG. 4. The same as Fig. 3, but forn > 1. The dotted, dashed, and solid lines are for m = 8, 10 
and 12 respectively. 

FIG. 5. The relation A and a for m = 6. In this case, n = 1. The dot-dashed line represents 
the COBE-DMR data at n = 1 and Qrms-PS = 15.3)uK, i.e., A ~ 3.5 x IQ-^. 

FIG. 6. The amplitudes of the power spectrum as a function of the thermal duration TV at 
three inflaton mass scales M and in the range of n < 1. The dotted, dotted and dashed lines are 
for m = 0, 2 and 4 respectively. 

FIG. 7. The same as Fig. 6, but for n > 1. The dotted, dotted and dashed lines are for 
m = 8, 10 and 12 respectively. 
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